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ABSTRACT 



P_|. We present IRTF SpeX observations of the M/L binary system 2MASS 

Q . J17072343-0558249. SpeX imaging resolves the system into a l'.'01±0.17 visual 

binary in which both components have red near infrared colors. Resolved low- 
■ resolution (R~150) 0.8-2.5 fxm spectroscopy reveals strong H2O, CO and FeH 

bands and alkali lines in the spectra of both components, characteristic of late- 
^ . type M and L dwarfs. A comparison to a sample of late-type field dwarf spectra 

indicates spectral types M9 and L3. Despite the small proper motion of the sys- 
tem (0'.'100±0'.'009 yr~^), imaging observations over 2.5 yr provide strong evidence 
that the two components share common proper motion. Physical association is 
also likely due to the small spatial volume occupied by the two components 
(based on spectrophotometric distances estimates of 15±1 pc) as compared to 
the relatively low spatial density of low mass field stars. The projected sepa- 
ration of the system is 15±3 AU, similar to other late-type M and L binaries. 
Assuming a system age of 0.5-5 Gyr, we estimate the masses of the binary com- 
ponents to be 0.072-0.083 and 0.064-0.077 Mq, with an orbital period of roughly 
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150-300 yr. While this is nominally too long a baseline for astrometric mass mea- 
surements, the proximity and relatively wide angular separation of the 2MASS 
J1707-0558AB pair makes it an ideal system for studying the M dwarf/L dwarf 
transition at a fixed age and metallicity. 

Subject headings: binaries: visual — stars: individual (2MASS J17072343-0558249) 
— stars: low mass, brown dwarfs 

1. Introduction 

Wide field photographic imaging and proper motion surveys of the past were largely 
incapable of producing an accurate census of nearby very low mass (VLM; M < 0.1 M©) 
stars and brown dwarfs. Late-type dwarfs emit weakly at visual bands, but are brighter at red 
optical and near infrared wavelengths. Progress in optical and infrared detector technology 
have paved the way for the current generation of wide field sky surveys, most notably the 
Deep Near Infrared Sky Survey (Epchtein et al. 1997), the Sloan Digital Sky Survey (York 
et al. 2000) and the Two Micron All Sky Survey (Skrutskie et al. 2006, hereafter 2MASS). 
Searches for faint red sources in these catalogs (e.g., Delfosse et al. 1997; Kirkpatrick et al. 
1999; Fan et al. 2000) have identified hundreds of very low mass stars and brown dwarfs, 
greatly expanding our VLM census and revealing the new spectral classes L (Kirkpatrick et 
al. 1999; Martm et al. 1999) and T (Burgasser et al. 2002; Geballe et al. 2002). 

The nascent field of VLM stars has been confronted with a number of fundamental 
questions regarding the properties of low luminosity sources, including formation scenar- 
ios, thermal evolution, chemical compositions and dynamics of cool stellar/substellar atmo- 
spheres, and the initial mass function. Multiple star systems are important laboratories for 
understanding these physical properties. Binary systems are key in the determination of 
stellar masses, which can be used to calibrate theoretical evolutionary and structure models. 
Binary parameters probe the star formation process, and relative comparisons can be made 
between coeval components. Searches for cool dwarf binaries through high resolution imag- 
ing have been conducted to explore the nature of VLM stars and brown dwarfs (e.g., Martin 
et al. 1999; Koerner et al. 1999; Reid et al. 2001; Close et al. 2002, 2003; Bouy et al. 2003; 
Burgasser et al. 2003a; Gizis et al. 2003; Siegler et al. 2005; Law et al. 2006; Burgasser et al. 
2006). The result of these efforts has been the discovery of roughly 75 VLM binary systems 
with varying mass ratios and projected separations (cf. Burgasser et al. 2006). 

The properties of the current sample of VLM stars and brown dwarfs suggests that both 
the multiplicity fraction and peak of the semimajor axis distribution are directly related to 
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the primary star mass. G to M binaries tend to peak at separations of 3-30 AU (Fischer & 
Marcy 1992; Henry, & McCarthy 1993; Reid & Gizis 1997), with binary fractions /^j^ ranging 
from ~65% (Duquennoy & Mayor 1991) down to ~30% (Fischer & Marcy 1992; Reid et al. 
2004; Delfosse et al. 2004). In contrast, nearly all VLM binaries have separations p < 20 
AU, with a binary fraction fun ~ 15% (e.g.. Close et al. 2003; however, see also Maxted & 
Jeffries 2005). The low frequency and preference of small separations for VLM binaries is a 
challenge for star formation theories, and a clear understanding of binary parameters as a 
function of mass, age, and metallicity provides empirical clues for understanding VLM star 
and brown dwarf formation processes. 

Studies of VLM binaries also reveal the detailed properties of cool star and brown 
dwarf atmospheres. Dust formation and evolution remains an outstanding problem across 
the spectral transition from M dwarfs to L dwarfs (hereafter the M/L transition), where 
lower temperatures and higher pressures enable some molecules to form solid condensates 
and descend from the photosphere, changing the overall morphology of the spectrum and 
the atmospheric pressure/temperature profile (Lunine et al. 1989; Tsuji et al. 1996; Lodders 
1999, 2002; Burrows & Sharp 1999; Ackerman & Marley 2001; Lodders & Fegley 2002). 
In particular, M dwarfs are characterized by their strong TiO and VO absorption bands; 
but as temperatures descend into the L dwarf regime, refractory elements such as Fe, Mg, 
Ti, V, Al, and Ca are removed from the gas in the photosphere by the condensation and 
sedimentation process. Theoretical models examining these processes are compared to the 
existing spectra of field dwarfs, which have a variety of metallicities, masses, and ages. M/L 
binaries, on the other hand, can help clarify the existence of atmospheric condensation by 
isolating parameters such as age and metallicity, assuming coevality. 

A total of 14 binaries comprised of M and L dwarf components have been discovered 
to date (Gizis et al. 2000, 2003; Bouy et al. 2003; Close et al. 2003; Freed, Close & Siegler 
2003; Martin et al. 2003; Siegler et al. 2003, 2005; Chauvin et al. 2004, 2005; Billeres et al. 
2005; Burgasser & McElwain 2006, see Burgasser et al. 2006). This paper presents the dis- 
covery of a new M/L binary\ 2MASS J17072343-0558249AB (hereafter 2MASS 1707-0558), 
identified via resolved near infrared imaging and spectroscopy using the SpeX instrument 
(Rayner et al. 2003) mounted at the 3 m NASA Infrared Telescope Facility (IRTF). In § 2 we 
describe our observations and data reduction, and present results. In § 3 we analyze these 
data, determining individual spectral types for the resolved components and the photometric 
properties of the system. We argue for physical association by common proper motion, sim- 
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ilar spectrophotometric distances, and comparing the volume occupied by the M/L pair to 
the measured density of VLM dwarfs in the Solar Neighborhood. We also make preliminary 
estimates of the individual masses and orbital characteristics of the system. We discuss our 
results in § 4, placing the 2MASS 1707-0558 system in context with other VLM binaries. 
This work is summarized in § 5. 

2. Observations 

The unresolved source 2MASS 1707-0558 was discovered by Gizis (2002) in a search for 
late-type dwarfs in the direction of the TW Hydrae association using the 2MASS survey. 
Optical spectroscopy of the composite system infers a spectral type of M9 on the Kirkpatrick 
et al. (1999) latc-M and L dwarf scale, and confirms this source as a normal field dwarf. Ha 
emission, common for late-type M dwarfs (Gizis et al. 2000; West et al. 2004), was observed 
with an equivalent width of 0.4 A. The signal-to-noise and resolution of the optical spectrum 
was not sufficient to detect the 6708 A Li I fine, so this source was suspected to be a nearby 
Hydrogen-burning star or brown dwarf. 

2.1. Imaging 

2MASS 1707-0558 was first observed on 2003 March 23 (UT) as a spectral comparison 
star for the 2MASS Wide Field T Dwarf Search program (Burgasser, McElwain & Kirk- 
patrick 2003b). Observing conditions were good with 0'/5 seeing at J-band, and the SpeX 
imager/guider was used to sample a 60"x60" field of view at 0'/12 pixel~^. While acquiring 
2MASS 1707-0558 with the imager, we resolved two sources at the given sky position. We 
subsequently obtained J, H and K"^ images of the visual pair for color comparison. In- 
tegrations of 15 s were obtained in 4 dithered exposures on the chip. A second series of 
images was obtained on 2004 August 9, with slightly hazy conditions and 0'.'9 seeing at J . 
Images were slightly out of focus during the second campaign. Integrations of 20, 15, and 
12 s with 4 dithers were obtained in the J, H, and iiT-band, respectively. We observed the 
single star USNO-A2.0 0825-10078125 (Monet et al. 1998) concurrently with the 2MASS 
1707-0558 observations to serve as a point spread function (PSF) calibrator. Finally, a third 
epoch of images was obtained on 2005 August 10, during poor conditions and 0'/9 seeing at 
J. Integrations of 10 s with 4 dithers were obtained in both the J and AT-band. In 2005, 



^JHK filters for SpeX are based on the Mauna Kea Observatories near-infrared (MKO-NIR) system 
(Simons & Tokunaga 2002; Tokunaga, Simons & Vacca 2002). 
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we observed the nearby star USNO-A2.0 0825-10079794 as a PSF calibrator. A log of our 
imaging observations is provided in Table 1. 

All imaging data was reduced in a typical procedure for near infrared imaging. Flat fields 
were constructed with sky frames, which were median combined, subtracted by a median 
combined dark frame of the same integration time, and normalized. An image mask was 
constructed to distinguish deviant pixels that were excessively bright in the dark frames and 
cold in the sky fiats. The target and PSF images were pair-wise subtracted, divided by the 
normalized fiat, cleaned by linear interpolation over the bad pixels, and added together by 
integer pixel shifts to match the peak flux position of the brighter source. 

Reduced J, H and i^-band images from the 2003 observations of 2MASS 1707-0558 arc 
displayed in Figure 1. The two sources are clearly resolved in all three bands and separated 
by roughly 1" on the sky. We discuss the relative fluxes and astrometry for the pair in § 3.1. 

2.2. Spectroscopy 

Near-infrared spectra of both components were acquired on 2003 May 23 (UT) with 
SpeX in prism mode, using the 0'.'5 slit. This conflguration yields a single-order, low- 
resolution (R~150) 0.8-2.5 jum spectrum with a dispersion of 20-30 A pixel"^ onto the 
Aladdin 3 1024x1024 InSb array. Conditions during the observations were good with seeing 
of 0'.'5 at J-band, and the target was observed at an airmass of 1.16. In order to obtain the 
spectrum of each component separately, we aligned the slit perpendicular to the orientation 
of the two sources, guiding on one component while observing the other through the slit. 
Based on the slit width, component separation and PSF full-widths at half-maximum, we 
estimate that only 4-6% of the light from the guiding component contaminated the slit. Ex- 
posures of 120/180s were taken as dithered pairs for the brighter/fainter component. The 
AO star HD 171149 was observed immediately after the target exposures at similar airmass 
(1.17), followed by internal flat-fleld and Ar arc lamps for pixel response and wavelength 
calibration. 

Data were reduced using the Spextool package version 3.2 (Gushing, Vacca, & Rayner 
2004). The raw target data were processed by performing linearity corrections, pair-wise 
subtraction, and division by a normalized flat field. The target spectra were then extracted 
using the Spextool default settings for point sources, and wavelength solutions were calcu- 
lated using the Ar arc calibration frames. Extracted spectra from the same source were 
scaled to match the highest S/N spectrum of the set, and the scaled spectra were median 
combined. The resultant spectrum was removed of telluric features, intrinsic AOV star lines. 
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and instrumental response signatures following the procedures of Vacca et al. (2003). A Vega 
model spectrum was employed to produce a true telluric spectrum from the AO standard. 
The model spectrum was shifted to match the radial velocity of the standard spectrum, and 
then the modified model spectrum was reddened and scaled to match the standard spec- 
trum. Line shape kernels were constructed using unresolved arc lines from the wavelength 
calibration observations. These kernels were then used to broaden the model spectrum line 
widths and smooth the model spectrum to match the observed resolution. The model spec- 
trum was completed by adjusting the H line strengths to reflect that of the observed AO 
standard. The corrected telluric spectrum is constructed by dividing the model spectrum by 
the AO spectrum. The corrected telluric spectrum is then multiplied by the target spectrum 
to produce the final fiux-calibrated spectrum. 

The reduced spectra of the two components are plotted in Figure 2. The spectrum of 
the bright component exhibits TiO absorption at 0.76, 0.82, and 0.84 //m, VO absorption 
at 1.05 yum, K I doublets at 1.17 and 1.25 fim, the onset of FcH absorption at 0.98, 1.19, 
and 1.58 /xm, CO at 2.3 jj,m, and strong H2O absorption at 1.4 and 1.9 /im. The spectrum 
of the faint component has diminished TiO and VO absorption features, but still shows the 
CrH, K I doublets, FeH, CO and H2O absorption features. In addition, the peak fiux of 
the spectral energy distribution for the faint component is also shifted redward to 1.28 /xm. 
These spectral features suggest that the bright component is a late M dwarf, consistent with 
the optical spectral type of Gizis (2002); while the faint component shows features which 
are indicative of an early L dwarf. We discuss the classification of these sources in further 
detail in § 3.2. 

3. Analysis 
3.1. PSF Fitting 

The PSF wings of the two sources at the position of 2MASS 1707-0558 are slightly 
blended in the 0'.'5, 0'.'9, and 0'.'9 seeing of the 2003, 2004, and 2005 image data, respectively. 
We therefore determined the relative fluxes and astrometry of the pair through a PSF fltting 
algorithm similar to that employed by Burgasser & McElwain (2006). For the 2004 and 
2005 images, we used a comparison star as the PSF model. For each filter, the components 
of the 2MASS 1707-0558 pair were first fit to a two-dimensional Gaussian to determine the 
approximate fiux center and a rough estimate of the peak amplitude. A synthetic image was 
created by shifting, scafing and adding the corresponding PSF star images to match each 
of the component positions and fiuxes. This model image was subtracted from the original, 
and the standard deviation of the residual image was used to measure the quality of the fit. 
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The model image was iteratively modified, changing the position of the primary, the position 
of the secondary, the peak flux of the primary and the peak flux of the secondary, in that 
order, to improve the fit until a minimum in the standard deviation of the residual image 
was achieved. This routine was performed on every pair-wise subtracted image for each filter 
in order to estimate the experimental uncertainty of the fits. For the 2003 data, which did 
not include a PSF comparison star, we used the same algorithm with a fixed width Gaussian 
profile to model the PSF. The astrometric results — measured separations (p) and position 
angle {(f)) - are listed in Table 2. Relative magnitudes and aggregate astrometric measure- 
ments are given in Table 3. The uncertainties in all of these measurements correspond to 
the 1(7 scatter in the values computed from all image frames. 

Relative magnitudes were measured using MKO-NIR filters, but composite systemic 
photometry is on the 2MASS system. In order to derive an apparent magnitude for each 
component, we converted the MKO-NIR relative magnitudes (AM^^^) to the 2MASS sys- 
tem as 

^j^2MA55 _ j^2MASS _ j^2MASS 

= Mf - Mj^^o + {M^^^^^ - Mf - (Mf^^^^ - (2) 
= AM^^o + - (3) 

where S = M^^^^^ — M^^*^ is the filter translation factor. These values were determined 
by calculating synthetic magnitudes from the spectra of the 2MASS 1707-0558 components, 
using MKO-NIR and 2MASS relative response curves (i?Af (A) = filter M transmission func- 
tion X optical response x telluric absorption at an airmass of 1) from S. Leggett (2004, 
private communication) and Cohen et al. (2003)^, respectively:^ 

M = -2.51ogio \ I I (4 

Here, /^^^" is the observed spectrum of the AOV star Vega from Hayes et al. (1985). The 
filter translation factors were measured to be Sj^^—Sj^a — —0.026, SH,b—SH,a — 0.013, and 
^K,b~^K,a — 0.002. Using these values, individual component magnitudes were computed 
from the composite 2MASS magnitudes for the 2004 and 2005 observations; mean values 
and their errors are listed in Table 3. 



^See Cutri et al. (2003), § IV. 4. a: http://www.ipac.caltech.edU/2mass/releases/allsky/doc/sec6_4a.html#rsr. 

^Note that we do not include an additional factor of A/hc in Eqn. 5 for the MKO-NIR filters to compensate 
for the photon counting properties of current detectors (this factor is already included in the response curves 
of Cohen et al. (2003)). This omission only leads to a ~0.01 mag offset in the derived magnitudes for 
late- type dwarf spectra (Stephens & Leggett 2004). We thank our referee for pointing out this detail. 
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3.2. Spectral Classification 

Spectral classification of stars is traditionally performed by comparison to a template 
of spectral standards, observed with the same dispersion and signal-to-noise ratio (Morgan 
et al. 1943). Unfortunately, a set of M and L dwarf spectral standards does not yet exist in 
the near infrared, although methods of near infrared classification have been considered by 
various groups (Tokunaga & Kobayashi 1999; Reid et al. 2001; Testi et al. 2001; Burgasser et 
al. 2002; Geballe et al. 2002; McLean et al. 2003). Here, we simply compare the component 
spectra of 2MASS 1707-0558 to previously observed SpeX data for the optical late M and L 
dwarf spectral standards (Burgasser et al. 2004; Cruz et al. 2006, in preparation). Figure 3 
overlays the normalized 2MASS 1707-0558 component spectra with those of the comparison 
stars to demonstrate the best spectral type fit for each component. The northern component 
spectrum matches the M9 optical spectral standard LHS 2924 (Probst & Liebert 1983; 
Kirkpatrick et al. 1991) quite well, with equivalent depths of the TiO and VO bands, and 
maintaining a similar shape over the entire spectrum. The southern component spectrum 
exhibits diminished TiO and VO features and increased H2O absorption, matching that of 
the L3 field dwarf 2MASS J08472872-1532372 (Cruz et al. 2003). While M and L spectral 
type classifications are generally defined at optical wavelengths, the excellent agreement 
between the 2MASS 1707-0558 near infrared spectra and those of the spectral comparators 
suggest that our classifications are consistent with the optical types and accurate to within 
0.5 subclass. 

An alternative method for near infrared classification is the calculation of spectral indices 
of diagnostic regions of the spectrum which are correlated with spectral types. We measured 
the Reid et al. (2001) H20'^ and H20^ indices, sampling the 1.3 and 1.5/im steam bands, 
respectively; and the Kl index (Tokunaga & Kobayashi 1999), a probe of H2O absorption at 
2.0//m. We then used the spectral type calibrations given in Reid et al. (2001), appropriate 
over optical spectral types M8 to L6, to derive classifications of M9 and L2 for the brighter 
and fainter components, respectively, with an uncertainty of 2 spectral types, as determined 
by the dispersions in each of the linear spectral type fits. Therefore, the spectral indices 
confirm the spectral types assigned by the spectral comparison method, but the large errors 
of this classification method prohibits a more precise estimate of the spectral types. 

3.3. Are 2MASS 1707-0558AB Gravitationally Bound? 

The angular proximity of the 2MASS 1707-0558 components as well as their similar 
brightnesses and spectral types, infer that the two components are gravitationally bound. 
One of the more conclusive tests of physical association is the detection of common space 
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motion. We examined the astrometry of the system over 3 epochs spanning 2.47 years to 
determine if the pair shared common proper motion. Table 2 hsts the astrometric results 
from the PSF fitting, as well as the expected astrometric results if the secondary object 
were an unrelated, unmoving background source. Figure 4 displays the expected separation 
and position angle of the system over time if the system is not gravitationally bound, with 
overlying data points for the measurements from the PSF fitting. A consistent separation and 
position angle over time confirms common proper motion. The proper motion of 2MASS 
1707-0558, p = 0'.'100±0'/008 at 88±10° (from the SuperCosmos Sky Survey; Hambly et 
al. (2001a,b,c)), is small and implies angular motion of only 0'/25 over our observational 
basehne, about 2 pixels at the SpeX plate scale. Our observations are not precise enough to 
establish common motion through the angular separation of the two components; however, 
our position angle measurements (accurate to 2-3°) can rule out a stationary background 
source at the 2.9a" (99.6% confidence) level. 

We compare our spectral types and photometric relative magnitudes derived by PSF 
fitting to the predicted relative magnitudes for an M9 and L3 binary system. Cruz et 
al. (2003) produces a polynomial Mj/spectral type relation (based on 2MASS photometry) 
using a sample of ultracool dwarfs (M6-L8) identified in the 2MASS catalog with independent 
parallax measurements. Assuming the stars are located at the same distance and of spectral 
type M9 and L3, the empirical polynomial fit gives AMj=1.20±0.35, consistent with our 
value of AMj=1.71±0.15. There is no parallax measurement for 2MASS 1707-0558, but 
spectrophotometric distances can be estimated by comparing the apparent magnitudes of 
each component to the absolute magnitudes typical for M9 and L3 dwarfs, as calculated 
using low mass absolute magnitude/spectral type relations (Dahn et al. 2002; Cruz et al. 
2003; Vrba et al. 2004). The two observed sources in this system lie at a mean distance of 
15 pc, with a 1 pc standard deviation in our measurements. 

But what if these two sources are nearby but unaffiliated late-type dwarfs? We also 
find that the chance alignment of two similarly classified low mass stars within 1" of each 
other is highly improbable. Cruz et al. (2003) performed a wide field search for VLM stars, 
discovering 30 M9-L3 stars over an area of 16,350 dcg^ at a depth of 20 pc, indicating a 
space density for these stars of p ~ 0.002 pc^'^. At a distance of 15±1 pc and an angular 
separation of 1", we estimate the volume of space occupied by the components of 2MASS 
1707-0558 toh&V 3.3x10"^ pc^. Assuming Poisson statistics, the probabihty (P) that 
one or more random low mass stars would be located within this volume of space can be 
calculated as 

P = 1 - e-^^, (5) 

or 10~^°, ruhng out chance ahgnment with high confidence. All the tests discussed in this sec- 
tion lead us to conclude that the components of 2MASS 1707-0558 comprise a gravitationally 
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bound system. 

4. Discussion 

4.1. 2MASS 1707-0558AB System Characteristics 

We calculate the projected physical separation of 2MASS 1707-0558AB to be 15±3 AU, 
based on the astrometric measurements and derived spectrophotometric distance. Temper- 
atures for the two components can be derived using the Te///spectral type relation for low 
mass stars from Golimowski et al. (2004), yielding Tg// = 2400±175 K for 2MASS 1707- 
0558A and 1950±190 K for 2MASS 1707-0558B. We calculate individual masses using the 
Burrows et al. (1997) models for ages of 0.5 and 5 Gyr, sampling the typical ages for field 
dwarfs (see below). The derived masses are 0.072-0.083 Mq and 0.064-0.077 Mq for the A 
and B components, respectively, with a total system mass of 0.136-0.160 M© and a mass 
ratio q = M2/M1 = 0.88-0.92. This mass ratio is again consistent with the properties of 
VLM binaries, which are predominantly near-equal mass systems. Combining the estimates 
of the physical separation and mass indicates an orbital period of roughly 150-300 yr, an 
unfortunately unreasonable timescale for dynamical mass measurements. A summary of the 
2MASS 1707-0558AB system characteristics can be found in Table 4. 

4.2. Broader Implications for M and L Dwcirfs 

This binary system is similar to other field binary systems, with a small separation 
and nearly equal mass components. 93% of known VLM binaries have separations less than 
20 AU, like 2MASS 1707-0558. and 77% have q > 0.8 (Burgasscr et al. 2006). Therefore, 
this system further strengthens the case that VLM binaries have specific traits that may be 
related to their formation mechanism. 

The M9/L3 binary 2MASS 1707-0558 is also a useful probe of the M/L transition, 
since the components have presumably coevolved with a fixed age and metallicity. Resolved 
binaries with spectroscopy are ideal systems to study the complex atmospheric chemistry 
of low mass stars, especially the thickening of dust, effects of condensation on the spectral 
energy distribution, and the relative heights of the cloud base. M dwarfs are characterized 
by their TiO and VO bands (Kirkpatrick et al. 1999), while the transition into the L dwarf 
regime is marked by the depletion of these molecules through condensation, leading to the 
expression of hydrides such as FeH and CrH (Lodders 2002). The atmospheric chemistry has 
been modeled for an assortment of field stars, yielding various masses, ages, and intrinsic 
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metallicities. A study of resolved binary systems straddling the M/L boundary, such as 
2MASS 1707-0558, will permit a fair comparison of low mass atmospheres. 

In addition, comparisons of the individual rotational velocities through atmospheric line 
broadening may constrain particular formation and evolution models. A direct measurement 
would determine if the components have equivalent or varied rotation rates, clarifying a 
dependency on either the mass or the origin of the object. Chromospheric activity is not 
well understood in low mass stars, but it is likely to be correlated to the spin rate of the 
star, possibly as a result of the magnetic fields induced by a turbulent dynamo mechanism 
(Durney et al. 1993, see also Chabrier & Kueker 2006). Field studies suggest that Ha activity 
increases from early-type to late-type M dwarfs, and then diminishes through the early-type 
L dwarfs (Gizis et al. 2000; Mohanty & Basri 2003; West et al. 2004). An unresolved optical 
spectrum of 2MASS 1707-0558 shows modest Ha emission (Gizis 2002), but the origins of 
this detection are unclear. Is this activity present in both components, and at what level 
does each component contribute to this detection? Furthermore, a measurement of the soft 
X-ray contributions of each component (resolvable with Chandra) will help further constrain 
the magnetic field mechanisms and their effects on coronal heating. 

The substellar nature of brown dwarfs can be directly tested through measurements of 
the 6708 A Li I line, but unfortunately the original optical spectrum was unable to resolve 
this line. The existence of primordial lithium is a conclusive determinant of substellar nature, 
as it is quickly destroyed when core temperatures reaches ~ 2x10^ K, a prerequisite reached 
by VLM stars with masses greater than ~0.06Mq (Rebolo et al. 1992; Magazzu et al. 1993). 
It should be noted that an object can lie below the hydrogen-burning minimum mass (i.e., 
substellar) and fail to exhibit Li 1 absorption. Depending on the age of the system, either one 
or both of the components will show Li 1 in absorption since they follow different evolutionary 
tracks. Hence, in the situation of a binary such as 2MASS 1707-0558 that may straddle the 
substellar boundary, the existence of the Li I line can be also used to trace the age of the 
system (cf., Liu & Leggett 2005). Three possible scenarios exist: 2MASS 1707-0558A and 
B could both be brown dwarfs and contain Li I undepleted from their atmospheres {t < 0.5 
Gyr), 2MASS 1707-0558A could be a star or brown dwarf with mass greater than ~O.O6M0 
and 2MASS 1707-0558B a Li-bearing brown dwarf with a mass less than ~O.O6M0 (0.5 Gyr 
< t < 1 Gyr), or both components could be low mass stars or brown dwarfs with masses 
greater than ^O.OGM© {t > 1 Gyr). 2MASS 1707-0558 is therefore a rare VLM field system 
that can be assigned an accurate age through the Li I diagnostic. 

An alternative age diagnostic, although less rehable, is a consideration of space kine- 
matics and chromospheric activity. Stars are generally born with the space motion of their 
natal cloud, and over time accumulate higher space velocities through interactions with the 
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components of the Galactic disk (Wielen 1977). Gizis et al. (2000) have found that a sample 
of field M and L dwarfs separate into two groups according to kinematics and chromosphcric 
activity. An old population of M dwarfs have high tangential velocities {ytan > 20 km s~^) 
and strong Hct emission, while a younger population exhibits low kinematics and low activ- 
ity. Assuming our spectrophotometric distance and the proper motion as measured by the 
SuperCosmos Sky Survey, 2MASS 1707-0558 has a Vtan of 7 km s~^. In comparison to the 
other late M dwarfs identified in Gizis et al. (2000), 2MASS 1707-0558 has one of the lowest 
tangential velocities and Ha measurements, suggesting that 2MASS 1707-0558 is associated 
with the young population. It is therefore likely to be closer to ~1 Gyr in age, but without 
additional constraints, we adopt a conservative age range of 0.5-5 Gyr. Figure 5 compares 
the properties of the 2MASS 1707-0558 components to evolutionary models from Burrows 
et al. (1997), clearly showing the secondary's possible location below the Li- burning mass 
limit. 



5. Summciry 

We have photometrically and spectroscopically resolved the 2MASS 1707-0558 into a 
M9/L3 binary separated by l'.'01±0.17. Physical association is confirmed by common proper 
motion, angular proximity, similar distances, and the statistical likelihood of two low mass 
stars occupying this small volume of space. System characteristics were derived by com- 
bining the imaging and spectral information obtained with the IRTF SpeX instrument, and 
considering the current evolutionary models and empirical large scale surveys. The derived 
J ,H jK-hand relative magnitudes are consistent with assigned spectral classifications and 
imply a spectrophotometric distance of 15±1 pc. The properties of 2MASS 1707-0558AB, 
including close separation and high mass ratio, are typical for late M/L binary systems. The 
angular separation of this system enables resolved spectroscopic measurements critical for 
the studies of low-mass star formation, atmospheric chemistry, and activity across the M/L 
transition. 

We thank our telescope operators Bill Golisch, Dave Griep, and Paul Sears, and instru- 
ment specialist John Rayner, for their support during the IRTF observations. We would 
also like to thank James Larkin, Stan Metchev, and Peter Plavchan for many useful conver- 
sations regarding the science presented herein, Michael Gushing for discussions on synthetic 
photometry, and our anonymous referee for a careful review of the manuscript. 
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Fig. 1. — The reduced mosaic discovery images of 2MASS 1707-0557AB in the J- (left), H- 
(middle), and X-bands (right) obtained on 2003 March 23. Images are 6" on a side, with 
North (arrow) and East (hne) indicated in the first panel 
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Fig. 2.— Near infrared SpeX spectra of 2MASS 1707-0558A (top) and B (bottom), with 
the major absorption features of TiO, VO, H2O, FeH, CrH, Na I, and K I indicated, as well 
as regions affected by telluric absorption (circled-plus symbols). All data are normalized at 
1.27/xm and offset by a constant. The zeropoint of each spectrum is designated by a dotted 
line. 
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Fig. 3 — Comparison of 2MASS 1707-0558A and B spectra (solid lines) to optically de- 
fined spectral standards (dashed lines). (Left) 2MASS 1707-0558A overlaid on the spectral 
standards VB 10 (M8; Probst & Liebert 1983; Kirkpatrick et al. 1991), LHS 2924 (M9; 
Probst & Liebert 1983; Kirkpatrick et al. 1991), and 2MASS 0345+2550 (LO; Kirkpatrick et 
al. 1997). (Right) 2MASS 1707-0558 overlaid on 2MASS 0847-1532 (L2; Cruz et al. 2003), 
SDSS 2028+0052 (L3; Hawley et al. 2002), and 2MASS 1104+1959 (L4; Cruz et al. 2003). 
All data are normalized at 1.27^m and offset by a constant. 



- 21 - 



1 00 n — ' — ' — ' — \ — ' — ' — ' — \ — ' — ' — ' — \ — ' — ' — ' — \ — ' — ' — ' — ] 1 5 F 




-200 D , , , \ , , , \ , , , \ , , , \ , , -10 L , , , \ , , , \ , , , \ , , , \ , , , 

200 40G 600 800 1 000 200 40G 600 800 1 000 

Time (days) Time (days) 



Fig. 4. — Astrometric motion of 2MASS 1707-0558B relative to A. Curves show the angular 
separation (left) and position angle (right) expected for an unmoving, unassociated back- 
ground source, with dotted lines corresponding to the la errors in the estimated positions. 
The straight line marks the position of a companion sharing common proper motion with the 
primary. Data points are the mean values from the PSF fitting algorithm, with the errors 
corresponding to the la uncertainty in these measurements. The periodicity of the relative 
background source position is due to the parallactic motion of the primary at its estimated 
spectrophotometric distance (15±1 pc). 
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Fig. 5. — T^eff evolution as computed by Burrows et al. (1997, solid lines) for masses of 10, 
20, 40, 50, 60, 70, 75, 80, and 100 Mj„p (1 Mjup ^ 0.001 Mq). The tracks are labeled from 
left to right and bottom to top, respectively. The mass/age loci of the 2MASS 1707-0558 A 
and B components are indicated, based on Te//S derived using the empirical Tg^j/spectral 
type relation of Gohmowski et al. (2004) and an age estimate of 0.5-5 Gyr. The L3 (B) 
component is hkely to be a brown dwarf, while the M9 (A) component is located near the 
hydrogen-burning minimum mass. The dotted line displays the boundary where more than 
90% of the primordial lithium abundance is depleted, according to the models of Burrows et 
al. (1997). 
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Table 1. Observing log. 



UT Date 


Observation 


tint (s) 


Airmass 


J-band Seeing (") 


Calibrator 


2003 March 23 


JHK imaging 


60, 60, 60 


1.17 


0.5 


none 




0.8-2. 5/im spectroscopy 


480 (A), 720 (B) 


1.16 


0.5 


HD 171149 


2004 August 9 


JHK imaging 


80, 60, 48 


1.17 


0.9 


USNO-A2.0 0825-10078125 


2004 August 10 


JK imaging 


40, 40 


1.15 


0.9 


USNO-A2.0 0825-10079794 



Table 2. 2MASS 1707-0558 System Astrometry. 



Measured Astrometry Expected Astrometry 
If Background Source 
Date Offset PA Offset PA 

(aJTcsec) (deg) (arcsec) (deg) 



2003 March 23 1.04±0.04 145±2 

2004 August 9 0.99±0.09 144±3 0.97±0.04 147±2 
2004 August 10 0.99±0.14 142±3 0.94±0.04 152±2 
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Table 3. 2MASS 1707-0558AB System Properties. 



Parameter 


Value 


Kei 


«»■ 




1 




-05° 58' 24^9 


1 




(y.'lOOiOyOOS yr-1 


2 


e 


88±10° 


2 




15±1 pc 


3 


p 


ivoiztC'iy 


3 




15±3 AU-^ 


3 


<t> 


145±3° 


3 


Jd 


12.05±0.02 mag 


1 


Hd 


11.26±0.03 mag 


1 


Kd 


10.71±0.02 mag 


1 


AJ 


1.71it0.15 mag 


3 


AH 


1.14±0.10 mag 


3 


AK 


1.18±0.12 mag 


3 


Mtot" 


0.136-0.160 Mq 


3,4 


(f 


0.88-0.92 


3,4 


Period*' 


~150-300 yr 


3,4 



^Equinox J2000 coordinates at epoch 1999 
April 09 from 2MASS. 

''Spectrophometric distances derived us- 
ing absolute magnitude/spectral type rela- 
tions; sec § 3.3. 

''Physical separation derived from the an- 
gular separation and the spectrophometric 

distance. 

d2MASS photometry of the unresolved 
system. 

''Assuming an age of 0.5-5 Gyr 

References. — (1) 2MASS (Skrutskie 
et al. 2006); (2) SupcrCosmosSky Survey 
(Hambly et al. 2001a,b,c); (3) This paper; 
(4) Burrows et al. (1997). 
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Table 4. 2MASS 1707-0558AB Individual Properties. 



Parameter 


2MASS 1707-0558A 
Value 


2MASS 1707-0558B 
Value 


SpT 






Te//"^ (K) 


2400±175 


1950±190 




12.26±0.11 


13.96±0.11 




11.59±0.07 


12.72±0.04 




11.03±0.08 


12.20±0.08 


Estimated Mass'^ (Mq) 


0.072-0.085 


0.064^0.077 



"Spectral type uncertainty of ±0.5 subclasses. 

''Derived from the Te///spectral type relation of Golimowski et al. 

(2004). 

"^Magnitude on the 2MASS photometric system. 

'^Estimated masses from 0.5 to 5 Gyr, based on the evolutionary models 
of Burrows et al. (1997). 



